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Abstract 

We present a lattice QCD calculation of the heavy quark expansion parameters /i 2 and [Aq for 
heavy- light mesons and heavy- light-light baryons. The calculation is carried out on a 20 3 x48 
lattice at /3 = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy 
quarks. We obtain the parameters ^ and [Iq in two different methods: a direct calculation of the 
matrix elements and an indirect calculation through the mass spectrum, and confirm that the both 
methods give consistent results. We also discuss an application to the lifetime ratios. 
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I. INTRODUCTION 

The heavy quark expansion (HQE) is a fundamental tool in the study of heavy 

quark physics. The inclusive decay rate of heavy hadrons containing a single heavy quark 
may be expanded in terms of inverse heavy quark mass I/tuq using the Operator Product 
Expansion (OPE) technique, which enables us to calculate the inclusive rates in a model 
independent manner U, l^l |(| . In particular, the determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements \V c b\ and \V u b\ through the corresponding semi-leptonic 
branching fractions relies on HQE. 

It requires, however, several nonperturbative parameters as coefficients in HQE. At the 
order 1/toq the nonperturbative parameters 
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£{H Q ) = — — (H Q QiiDfQ H Q ) , (1) 

appear in general. Here, Q denotes a heavy quark field defined in the Heavy Quark Effective 
Theory (HQET), and \Hq) represents a heavy-light meson or a heavy-light-light baryon 
state (for b hadrons, = B, B*, A&, etc.). Both parameters have mass dimension 

two, since they include a (spatial) covariant derivative squared D 2 or a chromomagnetic 
operator B. The inclusive decay rate of Hq is written in terms of (j, 2 (Hq) and ^q(Hq) as 



T(H Q - X 



1927T 3 



2^q / m Q 



(3) 



where the coefficients c{ and c{ are perturbatively calculable. On the other hand, the 
parameters ^(Hq) and ^Lq{Hq) have to be extracted from some experimental data or to be 
calculated nonperturbatively. Several methods to determine fi 2 and \Xq have been studied, 
and some of them are summarized in Section |Hj 

In this work we calculate /i^ and \Xq in quenched lattice QCD using the NRQCD action 
including 0(1/i71q) terms for heavy quark. Since the matrix element of power divergent 
operator Q(iD) 2 Q suffers from large perturbative uncertainty in the matching calculation 
with the continuum operator |7(, we consider their difference between different hadron states, 
like /x^(Ab) — fi 2 (B), in which the power divergence cancels. This kind of difference is also 
interesting in its own right, as it appears in the evaluation of lifetime difference of b hadrons 

0. 



One of the advantages of this calculation is that we can choose several quark masses in 
the calculation so that the heavy quark mass dependence of the hadron masses and matrix 
elements may be studied. We calculate both matrix elements /ij; and fi G and compare them 
with the corresponding mass spectrum and its heavy quark mass dependence. Another 
advantage in the use of the NRQCD lattice action is that the statistical signal in the Monte 
Carlo calculation is much better than in the static limit Q]. 

This paper is organized as follows. In Section |Hl the implications for the heavy quark 
expansion parameters from heavy hadron spectrum and the results of the previous nonper- 
turbative calculations are discussed. In Section IIHI we describe our lattice calculation in 
detail. The results for hadron masses and heavy quark expansion parameters are shown in 
Section I1V1 The consistency check between the calculation of matrix elements and spec- 
trum is also presented. Our results are applied to the lifetime ratio of different b hadrons in 
Section El The conclusions are given in Section IVII 



II. HEAVY QUARK EXPANSION PARAMETERS 



In this section we briefly review the determination of the HQE parameters from mass 
spectrum and from some nonperturbative techniques. The determination through the mea- 
surements of several mass and energy moments in the inclusive B — > XJv and B — > X s 7 



decays is another possibility 



flllfl, II, 12, JjJ 



which is not covered in the following. 



A. Implications from spectroscopy 

The HQE parameters /i^ and \j? G defined in (JTJ) and (j2J) can be indirectly obtained through 
heavy hadron masses, using the HQE of hadron masses 

2m Q \m 2 Q J 

where A is the residual energy difference between Mh q and tuq surviving in the infinite heavy 
quark mass limit. The parameters fi 2 and n G appear in the correction term of 0(l/mg). 
Considering proper mass differences, certain combinations of A, [v\ and /j, g can be extracted 
as shown below. 

The notation Ai and A2 is often used instead of [v\ and /j, g for B and B* mesons in the 
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literature. The relation between Ai )2 and /i 2 G is given by 
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and the HQE of meson masses in (jlj) becomes 



Ax = = -£(B*), (5) 

A 2 = l/4(B) = -vUB*), (6) 



M B =m b + A-±±lh + 0(±), (7) 
2rrih \ mi 1 



The parameter A 2 may be evaluated through the hyperfme splitting of ground state B mesons 
as 

M B * - M B ( ~ — ) = 46 MeV, (9) 
V 2m b J 

or, equivalently 

A 2 ~ i (m|» - M|) = 0.12 GeV 2 , (10) 

at the leading order. 

For A;, baryon, the parameter //^(A?,) vanishes, since the light degrees of freedom is spin 
singlet inside The relations 

= (ii) 

i/x 2 G (E 6 ) = -^(EJ), (12) 

hold for Sf, and baryons, as they are related by spin rotations, analogous to © and (jHJ) 
for B^ mesons. 

The spin-averaged meson mass becomes independent of A 2 

M g = -; = m b + A - - — + O — , 13 

but Ai cannot be extracted solely from this expression, as it appears together with the lowest 
order parameter A. In order to proceed further, we have to consider a similar relation for 
the D meson and take a mass difference to obtain 

M B - M B = m - m c - A, (-L - -L) + O (-U . (14) 
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The leading dependence on the heavy quark masses mj and m c can be subtracted out if we 
take a double mass difference 

M A b) - M^V - 2 i i H u 

M B M D 

from which we obtain 

Hl{k h ) - nl(B) = -0.01 ± 0.03 GeV 2 . (16) 

This argument relies on HQE truncated at order I/thq, which is questionable for charmed 
mesons and baryons. Therefore, for the use of the HQE parameter /i 2 in other phenomeno- 
logical analysis, some independent theoretical calculations are desirable. 




B. Nonperturbative calculations 

The determination of Ai using the QCD sum rule has been attempted by two groups 
and reached conflicting results A a = -0.5±0.2 GeV 2 [15J and -0.1±0.05 GeV 2 |lj. Their 
difference is explained to come from non-diagonal matrix elements like (B\Q(iD) 2 Q\B f ) , 
where B' is an excited state of B meson . Since there is no definite way to evaluate these 
matrix elements at present, it is not straightforward to improve the determination of Ai 
within the QCD sum rule technique. 

The lattice QCD can also be used to determine the HQE parameters. In the lattice 
calculation of the matrix element (B\Q(iD) 2 Q\B) the subtraction of quadratic divergence 
is essential, since otherwise the perturbative expansion to relate lattice and continuum op- 
erators poorly converges pj- First lattice calculation with such nonperturbative subtraction 
was done by Crisafulli et al. using the HQET on the lattice, which was updated in 
Gimenez et al. Q, and the result is Ai = 0.09 ± 0.14 GeV 2 . 

Another possible approach on the lattice is to fit the measured mass spectrum for various 
heavy quark masses with the mass relation (j3J) . Ali Khan et al. [n| performed such analysis 
for b flavored mesons and baryons using the lattice NRQCD for heavy quark. Their result 
is Ai = —0.1 ± 0.4 GeV 2 for B meson. Kronfeld and Simone performed similar analysis 
with a larger set of lattice data of heavy-light mesons, and quoted Ai = —0.45 ± 0.12 GeV 2 . 
The caption „ t * for b bary011 is available oel y from Ali Khan et a, Q. The y qU o t e d 
fi 2 n (A b ) = -1.7 ± 3.4 GeV 2 . 
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For the parameter A2, Gimenez et al. |l8| found Ao = 0.07 ± 0.01 GeV 2 from the direct 
calculation of the matrix element. Ali Khan et al. [l9j estimated \2(Bd) = 0.069 ± 0.019 
GeV 2 and \2{B S ) = 0.078 ± 0.012 GeV 2 from the hyperfine splitting measured on the lattice. 

The difference of A between several heavy hadrons is only estimated from the mass 
difference. Ali Khan et al. [l^ estimated A(A 6 ) - A(B) = 415 ± 156 MeV, A(E 6 ) - A(A 6 ) 
= 176 ± 152 MeV and A{B S ) - A(B d ) - 81 ± 31 MeV. 

In this work we calculate /i 2 and Hq on the lattice for ground state mesons and baryons. 
We use the both methods, namely the direct measurement of the matrix elements and the 
extraction from the heavy hadron spectrum. The difference of A is also evaluated from the 
mass difference. 



III. LATTICE CALCULATION 



In this section we present the details of our lattice calculation, which include the defini- 
tion of the NRQCD action, simulation parameters, and the method to extract the matrix 
elements. The matching of lattice operators onto their continuum counterpart is also dis- 
cussed. 



A. Lattice NRQCD 



We use the lattice NRQCD action [21 
action used in this work is the same as in 



22 for heavy quark. The particular form of the 
23, Q- 



£Wrqcd = J2 Q\ x ){$x,v - k q( x , y))Q(y)- 



(17) 



x.y 



The kernel to describe the time evolution of heavy quark is given by 



, 1 s^uivii- 



a5H s 



aH n \ n 



In J t 



(18) 



t+i \ / 1 

where the index to label the spatial coordinate is suppressed. The operator 8± is defined 

as 54 _) (x,?/) = 5 Xi -i m 8g t #, and 

A( 2 ) 



Hn EE 



277V 



8H EE -C B 



2m ( , 



-a 



B. 



(19) 
(20) 



A( 2 ) is a lattice covariant Laplacian 



A^Q(x) = £A 4 (2) Q( 



i=i 
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£ [^(x)Q(x + i) + t/J(a; - i)g(ac - i) - 2Q(x) 



i=i 



(21) 



and the chromo-magnetic field B is denned as the clover-leaf type on the lattice 



221. The 



parameter n in the evolution kernel (|18|) is a positive integer introduced to stabilize un- 
physical momentum modes j^,!^. With these definitions the lattice NRQCD action (|17|) 



deduces to the usual continuum NRQCD action 



pcont 
^NRQCD 



D, 



a-B 



Q 



(22) 



2.1/ " 2M 
in the limit of vanishing lattice spacing. 

The parameters appearing in the NRQCD action ()17j) . tuq and cb at this order, have 
to be matched onto their continuum counterparts using perturbation theory. The matching 
of heavy quark mass itlq is done through the calculation of hadron masses as described 
later. On the other hand, the one-loop calculation for eg is unfortunately not yet available, 
so we use the tree level value cb = 1. However, we apply the mean field improvement of 
the gauge link variable U^(x) — > U^(x)/uo [25j everywhere it appears, with m a mean link 
value defined through the plaquette expectation value uq = (gTrC/p). With the mean field 
improvement we expect that the tree level matching is reasonably good. Furthermore the 
final predictions for the matrix elements deduced from our analysis are given in the static 
limit, which is irrelevant to the parameter c#. 

The four-component heavy quark field h used to construct the hadron interpolating fields 
is related to the two-component nonrelativistic field Q through the Foldy-Wouthuysen-Tani 
(FWT) transformation 

' Q 
o 



h = R 



(23) 



V 



with the rotation matrix R given by 



R=l- 



7 • A 

2m 



at order I/ttlq. Our convention for the gamma matrices is 

" i r: ! I! • 7 — 





(24) 



(25) 



and the spatial covariant derivative is defined as 



AiQ(x) = l - [Ui(x)Q{x + 1) - - i)Q(x - i)] . (26) 



B. Simulation details 

Our calculation is carried out in quenched lattice QCD at (3 = 6.0 on a 20 3 x 48 lattice. 
Gauge configurations are generated with the single plaquette action, and 515 configurations 
are analyzed. 

The NRQCD action including 0(l/mg) described in the previous subsection is adapted 
for heavy quarks. Five heavy quark masses arriQ = 1.3, 2.1, 3.0, 5.0, and 10.0 are simulated 
to study the tuq dependence of hadron masses and matrix elements. The details on the 
parameters for heavy quark are shown in Table 



For light quarks, the 0(a)-improved Wilson action 26J] with the non-perturbatively tuned 
coefficient csw = 1-769 [27[ is used. Three hopping parameters K = 0.13331, 0.13384, and 
0.13432 are employed to extrapolate to the chiral limit K c = 0.135284(8). The inverse lattice 
spacing a -1 = 1.85(5) GeV is determined through the p meson mass m p = 770 MeV. 

The strange quark mass am s =0. 0460(22) is fixed using m^/m p =0.644 as an input. 

C. Hadron masses 

The hadron masses are measured through the asymptotic behavior of two-point functions 

C(J; t) = J2(J(x, t)J^\6, 0)) - e~ E -\ (27) 

X 

for sufficiently large time separation t. With the NRQCD action, for which the bare heavy 
quark mass is subtracted from the formulation, we obtain the binding energy E sim from 
the two-point function. The interpolating operator J is chosen such that it shares the same 
quantum number with the hadron of interest. The hadrons and their interpolating operators 
we consider in this work are the following. 

B = d lAl5 h, (28) 
B* = djih, (29) 
A b (s z = +1/2) = e abc (u a C l5 d b )h c p (30) 



A b (s z = -1/2) = t abc {u a C^d b )hl, (31) 
Y> b (s z = +1/2) = --L eabc (u a C l3 d b )hl + ^e abc (u a C^^d b ^j h c v (32) 

^. ( _ a „ll + 111 Jb \ , c , 1 , a ~*N 



S 6 (s 2 = -1/2) = -y -e Qfcc ^ a C " 2 ' V j ^ + ■^e abe (u a C l3 (r)hl, (33) 
E 6 *(s* = +3/2) = e abc (Vc^-^/) ^, (34) 



S*( Sz = +1/2) = Jjje^KC^)^ + ^e a 6c (u a C^^d b ) hi, (35) 



£*(s, = -1/2) = -±=e abc (Vtf^ipV) hi + ^e ab c{u a C lz d b )hl, (36) 
S*(s 2 = -3/2) = e abc (u a C^-^d b ) h c v (37) 



Although the notations motivated from the b hadron spectrum are used, we use them for 
general heavy quark mass we consider. The light quark fields u and d denote the relativistic 
up and down quark fields, respectively. The heavy quark field h has a subscript f or |, 
which represents its spin component in the z direction. We assume the Dirac representation 
of gamma matrices, and s z means the z component of the spin of baryons. The charge 
conjugation matrix C has a representation C = 7072- The superscript a, b or c denotes a 
color index of quarks. 

The smeared operator is used at the source in ()27j) to enhance the overlap with 
the ground state. It is defined such that the heavy quark field is smeared according to an 
exponential form e~ a ' rb around the light quark field fixed at the origin, r is a distance from 
the origin, and the parameters a and b are measured for the pion wave function. Thus, they 
depend on the light quark mass, as listed in Table HI Although it is not an optimal choice 
for heavy hadrons, the plateau is satisfactory as we demonstrate later. 

The hadron mass is obtained through the relation 

M had = (Z m M - E ) + E sim , (38) 

where Z m is the mass renormalization factor which relates the bare quark mass M with the 
pole mass and E is the energy shift of the heavy quark. These factors are perturbatively 
calculated at the one-loop level in [23] for our choice of heavy quark action. We summarize 
these factors in Table ITT1 
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D. Matrix elements 



To calculate the expansion parameters /j% and y? G from three-point functions, we construct 
a ratio 

R (J t t') = Zt,MZ,t)Oi(y,t)jW\d,0)) 



with Oi either the kinetic operator 



O w = Q(iD) 2 Q, (40) 



or the chromo-magnetic operator 



G = Q{B-B)Q. (41) 



The interpolating operator J is one of the operators listed in the previous subsection. The 
asymptotic behavior of the ratio yields the corresponding matrix element. We fix the position 
of the operator at t' = 9 and move the sink t. 

E. Operator renormalization 

The matching of the operators O n and G with their continuum counterpart is known 



28, 



only at the tree level except in the static limit, where one- loop coefficients are known 
The perturbative expansion is especially dangerous for the kinetic operator O n , since it mixes 
with lower dimensional operators QDqQ and QQ and thus power divergences appear. We, 
therefore, consider the differences of matrix elements with different hadron states, such as 
/x^(Afe) — fi^(B), in which the effect of mixed operators cancel at the leading order in I/tuq. 
The effect remains at finite values of 1/tuq, and hence we take the infinite heavy quark mass 
limit after measuring the differences at several values of rriQ. 

The other operator G does not mix with lower dimensional ones in the static limit. 
However, once the I/tuq correction is introduced, the mixing with and the other lower 
dimensional operators appears since the NRQCD action contains the a ■ B term. Hence, we 
again consider the difference among different hadron states to cancel the mixing contribution 
and take the infinite heavy quark mass limit. 

One-loop matching of the lattice operators to the corresponding continuum operators 

28j and the 



is known in the static limit. For the kinetic operator O n it is calculated in 
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multiplicative part is given as 

Z w = 1 + 0.0687 g 2 . (42) 

Since we apply the tadpole improvement using the plaquette expectation value and its effect 
is to multiply the link variable by 1/uq, the corresponding one-loop contribution -^g 2 has 
to be subtracted from the one- loop coefficient, and thus we obtain 

Z w = l- 0.0146 g 2 , (43) 

whose numerical value at (3 = 6.0 is 0.975 if we use the boosted coupling g 2 = g 2 /iiQ = 1.70. 
For the spin-chromomagnetic operator Oq the one-loop calculation is found in |29j as 

^ ■ .2 „2 



Z G = l + g 2 (-— In m 2 Q a 2 + 0.437j , (44) 

where itlq denotes the heavy quark mass arising from the continuum theory. The tadpole 
improvement amounts to multiply 1/uq and the one- loop coefficient is modified as 

,2 ( 3 i„ _2 2 



Z G = 1 + g 2 (- — In m 2 Q a 2 + 0.104 J , (45) 

and its numerical value is 1.12 for the b quark mass my, = 4.6 GeV. For both cases the 
tadpole improvement acts to greatly reduce the perturbative coefficients. 



IV. RESULTS 



In this section, we present the results for hadron masses and matrix elements. The heavy 
quark mass dependence of the matrix elements from the direct calculation is studied carefully 
by two methods. We also make a comparison between the results from the direct calculation 
and from the indirect calculation. All errors of measured quantities are estimated by the 
single elimination jackknife procedure. 



A. Hadron masses 

In Figures Q and 121 we show the typical effective mass plots for relevant mesons and 
baryons. The plateau is convincing for the B and B* mesons (Figure [IJ in the time region 
starting around t — 8, while it starts later in time for baryons (Figure |2J) and is dominated 
by statistical fluctuations after t = 20. We therefore fit the data in the time interval [10,20] 
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for mesons and in [12,20] for baryons. The results for the binding energy are summarized in 
Table EH 

Because the light quark mass dependence of the binding energy is well described by a 
linear function as shown in Figure E3 and El we can extrapolate (interpolate) the binding 
energy to the chiral limit (to the strange quark). The binding energy at the chiral limit and 
the strange quark is also presented in Table UTTl 

B. Matrix elements 

The ratio Ri(J; t, t') defined in (|39J) is shown as a function of t in Figure |S] for B and B* 
mesons. It shows a statistically cleanest data with heaviest light (K = 0.13331) and lightest 
heavy (aM = 1.3) quarks. The plateau is very convincing and appears earlier in Rq than 
in R n , and then we fit the data with a constant in the time interval [17,25] for R n or [14,25] 
for Rq- For other mass parameters the data are noisier, but we can identify the plateau in 
the same time interval. Similar plots for baryons (A&, £&, and ££) are shown in Figure El 
Since the statistical error dominates earlier in time we truncate the fit range at t = 23. The 
results for the matrix elements ^ and fi G are summarized in Table HVl and IVt respectively. 

From Figures ITHTOl we see that the light quark mass dependence of the matrix elements 
is mild though the statistical error grows as light quark mass decreases. We therefore take 
a simple linear fit in the light quark mass to obtain the results in the physical light quark 
mass. 

On the other hand, the heavy quark mass dependence of the matrix elements is significant 
as shown in Figures lTTHT^I In particular, the matrix elements fi G {B) and fi G (B*) in FigurelT2l 
are both positive at finite heavy quark masses, and hence do not respect the symmetry 
relation ~[/, G (B) = —n G (B*) given in (JBJ). This is due to the effects of operator mixing 
of Qa ■ BQ with spin singlet operators as mentioned in the previous section. The similar 
violation of the relation ()12j) is found in Figure El for the matrix elements of baryons 
/4(£ 6 ) and /4(££)- 

In order to extract the prediction in the static limit, where the symmetry relations have 
to be satisfied, we perform a fit of data in terms of a quadratic function in 1/Mg with a 
constraint known in the static limit. For mesons the constraint is (j5J) or (jUJ), while for baryons 
we may impose (fTT]) or (fT2*j) . The fitting curves describe the data well while satisfying the 
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constraints as shown in Figures ITTHT^l The bare matrix elements extrapolated to the static 
limit are listed in Tables IVT1 and fVTTl Since the chromomagnetic operator Qa ■ BQ does not 
receive the additive renormalization in the static limit, we may extract the physical result 
from these numbers. We obtain 

X 2 (B) (= = -/4(£*)) = 0.076(39) GeV 2 , (46) 

l&pb) = -2/4(££) = 0.23(11) GeV 2 , (47) 

after multiplying the renormalization factor Zq = 1.12 defined in (1451) . 

For the other matrix element /z 2 , the difference of the matrix elements between different 
heavy hadrons has to be considered in order to avoid the additive renormalization due to 
the mixing with lower dimensional operators. It also helps to reduce the statistical error as 
it correlates among different hadrons. The results are 

/i 2 (Afc) - ill (B) = -1.3(1.8) GeV 2 , (48) 
fil (E fe ) - fil (At) = -0.2(2.5) GeV 2 , (49) 

which include the multiplicative renormalization factor Z n = 0.975 as calculated in ()43j) . 
The SU(3) breaking /i 2 (£> s ) — ^1(B d ) has also a phenomenological importance, as it appears 
in the evaluation of the lifetime ratio T(B s )/r(B ( i). Our result is 

fil (B s ) - ^l(B d ) = 0.09(26) GeV 2 . (50) 

Another way to extract these physical quantities is to take the differences before extrap- 
olating the data to the static limit. As an example, we plot the difference of the matrix 
element /i 2 between At, baryon and B meson in Figure El Since each matrix element /i 2 (At) 
or /i 2 (-B) has a quite similar heavy quark mass dependence as seen in Figures HP and [T3l 
the heavy quark mass dependence almost cancels in the difference (Figure We fit the 
data with a linear function in 1/Mg and obtain 

/i 2 (At) - fil(B) = -0.01(52) GeV 2 , (51) 

in the static limit. This result is consistent with the previous analysis PH|) within one 
standard deviation. Since the heavy quark mass dependence is numerically better controlled 
in this method, we quote ()51|) as our final result, while taking the other to estimate systematic 
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uncertainty arising from the heavy quark extrapolation. The results for other differences of 

/4(£ 6 ) - fil(A b ) = 0.28(68) GeV 2 , (52) 
l£(B.) - nl(B d ) = 0.066(80) GeV 2 . (53) 

The same strategy - - differentiate then extrapolate — works even for fi G , since the 
additive renormalization at finite heavy quark masses mostly cancel in the differences like 
lx G {B*) - n%{B) or p%(J%) - n G (T, b ). Figure H^l shows the difference ~ Mg(^&) as a 

function of 1/Mg. We find that the heavy quark mass dependence is much milder than the 
individual matrix elements as shown in Figure^] This cancellation of the 1/Mg dependence 
is easily understood from Figure or El because the mass dependence is similar for all 
heavy hadrons. The results are 

\ 2 {B) = -\{ij? g {B*) - tf a {B)) = 0.094(19) GeV 2 , (54) 
/4(£ 6 ) = -\{nUn) - Vain)) = 0.147(60) GeV 2 , (55) 

which are consistent with the results obtained by taking the difference after the extrapola- 
tion, (jlHJ) and respectively. 

All these results are summarized in Table IVIII[ where "method 1" means our preferred 
method (differentiate-then-extrapolate) while "method 2" denotes the other (extrapolate- 
then-differentiate) . 

C. Heavy quark expansion parameters from mass differences 

The parameters A, /x 2 and fi G can also be indirectly obtained from hadron masses using 
the mass formula We use the hadron masses measured on the lattice to obtain the HQE 
parameters. 

We plot the mass difference M^ b — as a function of the spin-averaged meson mass 
inverse 1/Mg in Figure IT71 The change from rriQ to M§ is benign at this order because 
the difference between I/vtiq and 1/M§ is of order 1/rriQ which we neglect in this analysis. 
Fitting the data with a linear function of 1/Mg we obtain 

A(A 6 ) - A(B) = 428(68) MeV, (56) 
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from the intercept. This result is in good agreement with a previous lattice calculation by 
Ali Khan et ai, A(A 6 ) - A(B) = 415(156) MeV Q. Our result is slightly larger than the 
experimental value, which is about 310 MeV for bottom and charmed hadrons as plotted 
in Figure El by bursts. To draw a definite conclusion, however, we have to take account 
of several systematic errors. The finite volume effect is probably the most important one, 
because the physical extent of our lattice ~ 2 fm may not be large enough for baryons. 
The slope of the mass difference M^ b — yields 

*4(A 6 ) - fj*(B) = -0.38(47) GeV 2 , (57) 



the matrix elements (15 lj) and also with 
obtained from a combination {M\ b — 



which is compatible with the direct measurement o: 
the phenomenological estimate —0.01(3) GeV 2 ^ 
M B )-(M Ac -M B ). 

Similar analysis can be performed for M^ b — M\ b , which is plotted in Figure fTSl We 
obtain 

A(£ 6 )-A(A 6 ) = 96(96) MeV, (58) 
/i 2 (E 6 ) - /i 2 (A fe ) = 0.29(66) GeV 2 , (59) 

which are also consistent with the previous work A(£&) — A(A&) = 176(152) MeV and /i 2 — 
/4(A&) ~ Q. 

The strange-nonstrange mass difference M§ s —M% d is plotted in Figure IT91 It is interesting 
to see that the data agree well with the experimental value for and D( s ) mesons including 
the slope in 1/Mg. A linear fit gives 

A(B s )-A(B d ) = 90(7) MeV, (60) 
f4(B s ) - fil(B d ) = 0.056(42) GeV 2 , (61) 

which may be compared with A(B S ) - A(B d ) = 81(31) MeV and fil(B a ) - nl(B d ) = 
0.10(28) GeV 2 obtained in [lj. 

The hyperfine splitting in the mesons M#* — Afg d and Mb* — Mg s and in the baryons 
M£ — M-£ b is plotted in Figures I2TH I2T1 and 1221 respectively, as a function of l/M B . The 
numerical values at each quark masses are given in Table 11X1 where the statistical error in 
the hyperfine splittings is greatly reduced because it is highly correlated within the spin 
multiplets. 
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For the B — B* splitting ( Figure l2*Uj) we observe a linear behavior which is consistent with 
the expectation that the hyperfine splitting is proportional to 1/toq. The intercept at 1/Mb 
is, however, slightly negative. Since the hyperfine splitting is exactly zero in the static limit, 
we attempt a constrained fit with a linear and quadratic terms in 1/Mg, which is also shown 
in Figure 1201 It indicates that the quadratic term is not negligible and amounts about 5% 
at the B meson mass. From the coefficient of the linear term we obtain 

\ 2 (B d ) = 0.051(16) GeV 2 . (62) 

The similar analysis for B s gives 

X 2 (B S ) = 0.053(8) GeV 2 . (63) 

The data and fit curves are shown in Figure^J For the baryon hyperfine splitting M^* — M^ h 
shown in Figure [22j the statistical error is so large that the intercept of the linear fit is 
statistically consistent with zero. The slope yields 

/4(£*) - /4(£ 6 ) = -0.13(11) GeV 2 . (64) 

The experimental values of Mb* — Mb, Mb* — Mb, Mb* — Mb s and M^* — Ms c are also 
shown in Figures l20| I2T1 and l22l (Ms* — M^, b has not yet been measured.) The lattice data 
are significantly lower than these experimental results as in many other quenched lattice 
calculations. It is partly due to the fact that the spin-chromomagnetic interaction term in 
the lattice NRQCD action ()20j) is matched to the continuum full theory only at the tree 
level, although the mean field improvement is applied. Another important uncertainty is in 
the quenching approximation, whose effect is not yet entirely uncovered. 

The numerical results given in this subsection are also summarized in Table lyTTTl together 
with the results from other groups and the experimental values. 

D. Consistency among matrix elements and mass differences 

Results presented so far indicate that the HQE parameters are determined consistently 
with the direct measurement of the matrix elements and with the indirect measurement 
through the mass differences. However, more stringent test is possible using the data at 
fixed light quark mass, whose statistical error is smaller than in the chiral limit. Although 
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the numerical values are unphysical, there is nothing wrong in the consistency check. For 
this purpose we use the data at K = 0.13331, which corresponds to the heaviest light quark 
mass. 

From (0) and (|SJ) the hyperfine splitting Mb* — Mb is given by 4:X 2 /2mb, or up to higher 
order 1/mj corrections, 

M 2 B *-Ml = -&{B*)+^ G {B). (65) 

In Figure ESI we plot the results for —^(B*) + Hq(B) as a function of l/M B together with 
the lattice measurement of M B „ — M B . We observe that the relation (J55j) is satisfied well 
in the heavy quark mass region 1/Mg < 0.2 GeV -1 . Towards lighter heavy quark mass the 
data deviate from the relation (|65J1 . which is an indication of higher order effect. Similar 
analysis can be done for the hyperfine splitting of heavy-light-light baryon, i. e. the E£ — 
splitting. Figure |2*%I shows the mass difference and the matrix element —A/j,q. Both are in 
good agreement within the large statistical error in the hadron mass measurement. 
The heavy-light meson-baryon mass difference M\ b — M B is given as 

M Ab -M B = A(A 6 ) - A(B) + -L U(A b ) - ^(B)] . (66) 

In Figure EI3 we plot M\ b — M B as a function of 1/Mg. The slope obtained from the fit of 
the mass difference yields an indirect estimate of /x^(A&) — fi*(B) as —0.03 ±0.15 GeV 2 . Our 
results for the direct measurement of /x^(A&) — fJ%(B) are plotted in Figure ESI where the 
indirect measurement is shown by a band. Both measurements are completely consistent 
with each other. 



V. LIFETIME RATIO: A PHENOMENOLOGICAL APPLICATION 



In the ratio of lifetimes of different b hadrons and the hadronic matrix elements 



\i v and jj,Q appear as 



r(H^) , , nl(HV) - fi{Hj») + cg /40gg) - f*UHj 2) ) +0 (—) 7 (67) 



mi \ ml 



with a perturbative coefficient Cq — 12 Q]. Our calculation of the differences of the matrix 
elements ^{H^) — ) and imq(H^) — iiq(H^ ) may be directly used to evaluate the 

lifetime ratios at the order \ jm\. 
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Using our results /4(A 6 ) -fil(B d ) = -0.01(52) GeV 2 and /4(A 6 ) ~n%{B d ) (= -3\ 2 (B d )) 
= —0.282(59) GeV 2 , which are from the direct calculation (method 1), the lifetime ratio of 
A b and B d is evaluated as 

444 = 0-984 ±0.012 ±0.003 + [ -M , (68) 

with mj = 4.6 GeV, where the first and second error comes from the statistical error of 
u^YAfe) — fi^(B d ) and /^(Ab) — //^(S^), respectively. As discussed in the previous works 
|l4| it may not explain the experimental value 0.76(5) unless the higher order effect in 
the l/m& expansion has a substantially large effect. Our calculation does not imply such 
a large correction to the matrix element /x 2 as shown in Figure I2H1 At the order 1/mf 
the spectator effect arises, for which the hadronic matrix elements of higher dimensional 



operators are necessary 



14j. A lattice calculation |30j of those matrix elements suggests 



that the spectator effects are indeed significant but do not appear to be sufficiently large to 
account for the full discrepancy. 

The lifetime ratio of B s and B d is obtained as 

4^4 = 1-001 ± 0.002 ± 0.002 + O (Ar ) , (69) 
r(B d ) \ml)> 1 ' 

using our results for /4(5 S ) - /i 2 (5 d ) =0.066(80) GeV 2 and /4(£ s ) - 
fj,Q(B d )=— 0.012(32) GeV 2 . This result may be compared with the experimental value 
r(B s )/r(B d ) = 0.949±0.038 



VI. CONCLUSIONS 



In this article we present a lattice QCD calculation of the heavy quark expansion param- 
eters A, /i 2 and \0q for the heavy-light mesons and heavy-light-light baryons. The lattice 
NRQCD action is used for heavy quark and the results in the static limit are obtained by 
an extrapolation. 

For /i 2 and fig, we performed a direct calculation of the matrix elements through the 
three-point functions. While the light quark mass dependence of the matrix elements is 
small, the heavy quark mass dependence is significant due to the effect of the additive 
renormalization. The large heavy quark mass dependence mostly cancels by considering 
the difference of the matrix elements between different heavy hadron states, in which the 
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additive renormalization cancels. We also estimate the differences of the HQE parameters 
by studying the mass differences between several heavy hadrons. 

We find that the lattice measurements of the matrix elements and n 2 G are consistent 
with the mass relations predicted by the heavy quark expansion. Our numerical results 
for the defferences of fj% in the heavy quark mass limit are compatible with the previous 
determinations from the meson mass spectrum. The deficit of the hyperfine splitting — the 
well-known problem of the quenched lattice calculation — is also reproduced in the direct 
calculation of the matrix element /Iq. 

A direct phenomenological application of our results is the evaluation of the lifetime 
ratios at the order l/raj). Previously such analysis implicitly assumed that the heavy quark 
expansion truncated at I /ml is valid down to the charm quark mass, as the parameter 
was determined using the combined mass difference including charmed mesons and baryons. 
Through the direct lattice calculation we have confirmed that such analysis is justified. The 
problem of the small lifetime ratio r{k h )jr{B^) still remains. 
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ttMn 



n 
a 
b 



1.3 



2.1 



3.0 5.0 10.0 



3 3 2 2 2 
0.2248 0.2530 0.2711 0.3074 0.3425 
1.2484 1.1840 1.1465 1.0794 1.0294 



TABLE I: Simulation parameters. The parameter a and b is for the smeared source e a ' r . 



aM 


n 


A B aE Z m aA 


1.3 


3 


0.547 0.914 0.140 1.234 1.464 


2.1 


3 


0.754 0.578 0.193 1.148 2.218 


3.0 


2 


0.855 0.381 0.219 1.097 3.072 


5.0 


2 


0.946 0.176 0.242 1.045 4.983 


10.0 


2 


1.011 0.040 0.259 1.010 9.841 



TABLE II: Perturbative factors to obtain the hadron mass from E s i m using (|38|) . The perturbative 
expansions are given as ciEq = a s A and Z m = 1 + a s B where the coefficients A and B are given 



m 



23j. For the numerical analysis we use a renormalized coupling ay(l/a) = 0.256 for a s at (5 
6.0. aA in the last column is defined as aA = Z m aMo — qEq. 
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aM 


K 


aE sim (B) 




aE sim (A b ) aE sim (E b ) 


dEsimi^b) 


1.3 


0.13331 


0.4928(15) 


0.5156(17) 


0.8101(54) 0.8573(60) 0.8631(63) 


2.1 




0.5145(17) 


0.5298(19) 


0.8256(60) 0.8725(65) 0.8766(67) 


3.0 




0.5247(19) 


0.5357(21) 


0.8324(66) 0.8787(70) 0.8816(72) 


5.0 




0.5327(24) 


0.5391(25) 


0.8386(92) 0.8830(84) 0.8843(86) 


10.0 




0.5376(37) 


0.5401(38) 


0.847(16) 0.891(14) 


0.891(14) 


1.3 


0.13384 


0.4754(18) 


0.4987(20) 


0.7680(73) 0.8210(85) 0.8275(91) 


2.1 




0.4976(21) 


0.5132(22) 


0.7849(83) 0.8366(92) 0.8411(97) 


3.0 




0.5083(23) 


0.5194(25) 


0.7927(95) 0.843(10) 


0.847(10) 


5.0 




0.5166(29) 


0.5229(30) 


0.801(12) 0.850(12) 


0.852(12) 


10.0 




0.5218(43) 


0.5241(44) 


0.813(21) 0.862(19) 


0.862(19) 


1.3 


0.13432 


0.4599(24) 


0.4836(26) 


0.728(12) 0.786(15) 


0.793(15) 


2.1 




0.4825(27) 


0.4983(29) 


0.746(14) 0.801(16) 


0.806(16) 


3.0 




0.4934(30) 


0.5047(32) 


0.756(15) 0.808(17) 


0.812(17) 


5.0 




0.5021(36) 


0.5085(38) 


0.767(19) 0.820(20) 


0.821(20) 


10.0 




0.5079(57) 


0.5099(55) 


0.782(32) 0.831(24) 


0.830(23) 


1.3 


K s 


0.4826(17) 


0.5057(19) 


0.7850(66) 0.8357(78) 0.8419(82) 


2.1 




0.5046(20) 


0.5201(21) 


0.8015(75) 0.8510(83) 0.8553(87) 


3.0 




0.5151(22) 


0.5262(23) 


0.8089(84) 0.8575(90) 0.8606(93) 


5.0 




0.5233(27) 


0.5296(28) 


0.816(11) 0.864(11) 


0.865(11) 


10.0 




0.5284(41) 


0.5307(42) 


0.827(18) 0.873(16) 


0.873(16) 


1.3 


K c 


0.4290(33) 


0.4534(36) 


0.652(17) 0.721(21) 


0.729(22) 


2.1 




0.4524(37) 


0.4688(39) 


0.673(20) 0.736(23) 


0.742(24) 


3.0 




0.4639(42) 


0.4756(43) 


0.684(23) 0.744(24) 


0.749(25) 


5.0 




0.4733(48) 


0.4796(52) 


0.699(28) 0.761(29) 


0.762(29) 


10.0 




0.4799(77) 


0.4814(74) 


0.722(49) 0.777(36) 


0.776(34) 



TABLE III: Binding energy of heavy-light mesons and heavy-light-light baryons. 
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aM 


K 


a 2 /4(#) 


a 2 »l(B*) 


a 2 *4(A*) 


a 2 d^ b ) 


a 2 ^* b ) 


1.3 


0.13331 


-0.2507(24) 


-0.2643(25) 


-0.2503(62) 


-0.2416(56) 


-0.2449(55) 


2.1 




-0.0998(51) 


-0.1075(53) 


-0.096(15) 


-0.083(11) 


-0.085(11) 


3.0 




-0.0546(90) 


-0.0558(94) 


-0.051(27) 


-0.036(19) 


-0.037(19) 


5.0 




-0.016(23) 


-0.003(24) 


-0.072(63) 


-0.053(52) 


-0.048(52) 


10.0 




-0.030(88) 


+0.018(88) 


-0.28(27) 


-0.35(25) 


-0.34(25) 


1.3 


0.13384 


-0.2525(30) 


-0.2663(30) 


-0.2545(81) 


-0.2416(82) 


-0.2427(82) 


2.1 




-0.1027(63) 


-0.1104(66) 


-0.101(19) 


-0.081(17) 


-0.081(16) 


3.0 




-0.059(11) 


-0.059(12) 


-0.058(36) 


-0.035(28) 


-0.036(28) 


5.0 




-0.022(28) 


-0.003(29) 


-0.103(89) 


-0.082(77) 


-0.076(79) 


10.0 




-0.05(11) 


+0.01(11) 


-0.39(39) 


-0.51(38) 


-0.52(39) 


1.3 


0.13432 


-0.2537(40) 


-0.2678(40) 


-0.257(13) 


-0.229(14) 


-0.225(14) 


2.1 




-0.1042(84) 


-0.1128(88) 


-0.099(31) 


-0.057(27) 


-0.054(26) 


3.0 




-0.062(15) 


-0.062(16) 


-0.049(58) 


-0.008(47) 


-0.005(47) 


5.0 




-0.029(37) 


-0.004(38) 


-0.13(15) 


-0.10(13) 


-0.09(14) 


10.0 




-0.09(14) 


+0.00(14) 


-0.45(66) 


-0.71(69) 


-0.76(70) 


1.3 


K s 


-0.2517(28) 


-0.2654(28) 


-0.2526(74) 


-0.2399(73) 


-0.2415(72) 


2.1 




-0.1014(58) 


-0.1092(61) 


-0.098(18) 


-0.078(14) 


-0.079(14) 


3.0 




-0.057(10) 


-0.058(11) 


-0.053(33) 


-0.032(25) 


-0.032(25) 


5.0 




-0.020(26) 


-0.003(27) 


-0.090(80) 


-0.069(68) 


-0.063(70) 


10.0 




-0.05(10) 


+0.01(10) 


-0.34(34) 


-0.45(34) 


-0.46(34) 


1.3 


K c 


-0.2567(55) 


-0.2712(55) 


-0.265(19) 


-0.227(20) 


-0.218(20) 


2.1 




-0.109(12) 


-0.118(12) 


-0.107(44) 


-0.049(40) 


-0.041(40) 


3.0 




-0.069(21) 


-0.067(21) 


-0.057(82) 


-0.001(70) 


+0.004(71) 


5.0 




-0.040(51) 


-0.006(52) 


-0.19(22) 


-0.15(20) 


-0.14(21) 


10.0 




-0.14(20) 


-0.01(20) 


-0.64(96) 


-1.01(98) 


-1.1(1.0) 



TABLE IV: Matrix elements fi 2 , for heavy-light mesons and heavy-light-light baryons. 
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aM 


K 




a 2 /4(5*) 


a 2 fi 2 G {A b ) 


a 2 /4(£ 6 ) 


a 2 /4(££) 


1.3 


0.13331 


0.2507(20) 


0.1840(10) 


0.2027(12) 


0.2190(22) 


0.1982(13) 


2.1 




0.2340(28) 


0.1588(12) 


0.1798(12) 


0.1979(28) 


0.1736(15) 


3.0 




0.2067(35) 


0.1258(15) 


0.1482(13) 


0.1674(34) 


0.1402(18) 


5.0 




0.1671(56) 


0.0827(21) 


0.1056(17) 


0.1273(48) 


0.0944(25) 


10.0 




0.112(11) 


0.0395(39) 


0.0581(23) 


0.0856(86) 


0.0427(44) 


1.3 


0.13384 


0.2515(24) 


0.1837(11) 


0.2033(15) 


0.2208(29) 


0.1976(16) 


2.1 




0.2350(33) 


0.1585(15) 


0.1804(15) 


0.2000(36) 


0.1732(19) 


3.0 




0.2079(43) 


0.1253(18) 


0.1488(17) 


0.1700(43) 


0.1397(23) 


5.0 




0.1680(68) 


0.0823(26) 


0.1060(22) 


0.1303(62) 


0.0936(32) 


10.0 




0.111(13) 


0.0399(46) 


0.0584(29) 


0.090(11) 


0.0407(55) 


1.3 


0.13432 


0.2526(29) 


0.1835(15) 


0.2031(20) 


0.2240(43) 


0.1970(23) 


2.1 




0.2363(42) 


0.1581(19) 


0.1806(23) 


0.2040(53) 


0.1729(28) 


3.0 




0.2094(55) 


0.1249(23) 


0.1488(25) 


0.1744(63) 


0.1395(34) 


5.0 




0.1694(88) 


0.0820(33) 


0.1057(32) 


0.1346(91) 


0.0931(46) 


10.0 




0.110(17) 


0.0404(59) 


0.0589(43) 


0.095(15) 


0.0389(77) 


1.3 


K s 


0.2512(22) 


0.1839(11) 


0.2030(14) 


0.2203(27) 


0.1978(15) 


2.1 




0.2346(31) 


0.1586(14) 


0.1801(14) 


0.1994(33) 


0.1734(18) 


3.0 




0.2075(40) 


0.1255(17) 


0.1484(15) 


0.1692(40) 


0.1399(21) 


5.0 




0.1677(64) 


0.0825(24) 


0.1058(21) 


0.1293(57) 


0.0939(29) 


10.0 




0.112(13) 


0.0398(44) 


0.0583(27) 


0.088(10) 


0.0415(50) 


1.3 


K c 


0.2544(40) 


0.1830(20) 


0.2037(29) 


0.2277(60) 


0.1959(33) 


2.1 




0.2384(56) 


0.1575(25) 


0.1816(33) 


0.2087(75) 


0.1722(40) 


3.0 




0.2117(75) 


0.1242(31) 


0.1497(36) 


0.1797(91) 


0.1388(48) 


5.0 




0.171(12) 


0.0813(44) 


0.1061(46) 


0.141(13) 


0.0919(66) 


10.0 




0.108(23) 


0.0412(78) 


0.0594(60) 


0.102(22) 


0.035(11) 



TABLE V: Matrix elements fifj for heavy-light mesons and heavy-light-light baryons. 
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aM 


K 


a 2 nl(B) = a 2 /4(£*) 


a 2 nl{A b ) 


a 2 /4(£&) = a 2 /4(EJ) 


static 


0.13331 


-0.057(61) 


-0.25(17) 


-0.22(13) 


static 


K s 


-0.061(69) 


-0.29(21) 


-0.26(17) 


static 


K c 


-0.09(13) 


-0.48(53) 


-0.54(50) 



TABLE VI: Matrix elements /x 2 for heavy-light mesons and heavy-light-light baryons in the static 
heavy quark limit. 



aM 


K 


a^h(B) a^UB*) a^l{A b ) a 2 ^(S 6 ) a^PS) 


static 


0.13331 


0.065(15) -0.022(05) 0.039(11) -0.020(06) 


static 


K s 


0.064(17) -0.021(06) 0.043(13) -0.021(07) 


static 


K c 


0.060(30) -0.020(10) 0.059(29) -0.029(15) 



TABLE VII: Matrix elements for heavy-light mesons and heavy-light-light baryons in the static 
heavy quark limit. 





Direct, calc. 
(method 1,2) 


Mass difference 


Other works 


Exp. 


A(A fe ) - A(B) [MeV] 
A(E 6 ) - A(A 6 ) [MeV] 
A(B S ) - A{B d ) [MeV] 




428(68) 
96(96) 
90(7) 


415(156) [19] 
176(152) [19] 
81(31) [12] 




A(A b )-ti(B) [GeV 2 ] 
^ b )-ti(A b ) [GeV 2 ] 
ti{B s )-ti{B d ) [GeV 2 ] 


-0.01(52), -1.3(1.8) 
0.28(68), -0.2(2.5) 
0.066(80), 0.09(26) 


-0.38(47) 
0.29(66) 
0.056(42) 


[19J 

o [is] 

0.09(4) [18], 0.10(28) [19] 


-0.01(3) 
0.06(2) 


\ 2 {B d ) [GeV 2 ] 
\ 2 {B S ) [GeV 2 ] 
/£(E 6 ) [GeV 2 ] 


0.094(19), 0.076(39) 
0.090(10), 0.082(22) 
0.147(60), 0.23(11) 


0.051(16) 
0.053(8) 
0.09(7) 


0.070(15) [18], 0.069(19) [19] 
0.078(12) [19] 


0.12(1) 



TABLE VIII: Results for the HQE paremeters. 
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aM 


K 


a{M B * - M B ) 


a(M s * -M s J 


a(M Ab - M B ) 


a(M St - M Ab ) 


1.3 


K s 


0.0231(08) 


0.0062(21) 


0.2851(63) 


0.055(09) 


2.1 




0.0155(07) 


0.0043(18) 


0.2852(72) 


0.052(11) 


3.0 




0.0111(07) 


0.0031(17) 


0.2855(81) 


0.051(12) 


5.0 




0.0063(07) 


0.0014(18) 


0.288(11) 


0.048(15) 


10.0 




0.0024(08) 


-0.0003(27) 


0.297(18) 


0.046(24) 


1.3 


K c 


0.0244(16) 


0.0081(61) 


0.205(17) 


0.075(28) 


2.1 




0.0164(14) 


0.0058(53) 


0.209(20) 


0.067(31) 


3.0 




0.0116(14) 


0.0041(50) 


0.212(23) 


0.063(34) 


5.0 




0.0063(15) 


0.0018(53) 


0.221(28) 


0.063(41) 


10.0 




0.0016(16) 


-0.001(20) 


0.241(49) 


0.055(61) 



TABLE IX: Mass difference between heavy hadrons. 
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FIG. 1: Effective mass plot for the B (top panel) and B* (bottom panel) mesons at K=0. 13331 
and aMo=1.3. Solid lines represent the fitting result with an error band of one standard deviation. 
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FIG. 2: Effective mass plot for the (top panel) and £& (middle panel) and (bottom panel) 
baryons at K=0. 13331 and aMo=1.3. Solid lines represent the fitting result with an error band of 
one standard deviation. 
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FIG. 3: Binding energy of the B and B* mesons as a function of light quark mass at aMo=1.3 
(top panel) and 10.0 (bottom panel). 
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FIG. 4: Binding energy for the and baryons as a function of light quark mass at aMo=1.3 

(top panel) and 10.0 (bottom panel). 
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FIG. 5: Ratio Ri(J;t,t' = 9) for (top panel) and for pr G (bottom panel) at K=0. 13331 and 
aM=1.3. Open (filled) symbols are the data for the B (B*) meson. Solid lines represent a constant 
fit with an fit interval [17,25] for or [14,25] for p? G . 
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FIG. 6: Ratio Ri(J;t,t' = 9) for //^ (top panel) and for /j,q (bottom panel) at K=0. 13331 and 
aM=1.3. Open circles, open triangles and filled triangles are data for A;,, T,(, and baryons, 
respectively. Solid lines represent a constant fit with an fit interval [17,23] for [j% or [13,23] for /j,q. 
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FIG. 7: Matrix element //^ for the B and B* mesons as a function of light quark mass at aMo=1.3 
(top panel) and 10.0 (bottom panel). 
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FIG. 8: Matrix element /j/q for the B and B* mesons as a function of light quark mass at aMo=1.3 
(top panel) and 10.0 (bottom panel). 
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FIG. 9: Matrix element ^ for the £5 and baryons as a function of light quark mass at 
aMo=1.3 (top panel) and 10.0 (bottom panel). 
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FIG. 10: Matrix element p? G for the A&, and baryons as a function of light quark mass at 
aMo=1.3 (top panel) and 10.0 (bottom panel). 
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FIG. 11: Matrix element ^ for the B and B* mesons as a function of l/Mg. The value in the 
static limit is obtained from a fit in terms of a quadratic function in 1/Mg with the constarait 
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FIG. 12: Matrix element [1q for the B and B* mesons as a function of 1/Mg. The values in the 
static limit are obtained from a fit in terms of a quadratic function in \/Mg with the constarait 
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FIG. 13: Matrix element ^ for the A&, and baryons as a function of 1/Mg. 
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FIG. 15: Difference of the matrix element /^(A^) — ^(B^) as a function of 1/Mg. Open circle 
denotes the result from method 2. 
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FIG. 16: Difference of the matrix element ^(X^) — ^(Sft) as a function of 1/Mg. Open circle 
denotes the result from method 2. 
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FIG. 17: M\ b — Mb as a function of l/Mg. The light quark mass is exptapolated to the chiral 
limit. 
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FIG. 18: Mg b — Ma 6 as a function of l/Mg. The light quark mass is exptapolated to the chiral 
limit. 
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FIG. 22: Hyperfine splitting Ms* — Ms b as a function of 1/Mg. The light quark mass is exptap- 



olated to the chiral limit. 
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FIG. 23: Hyperfine splitting of the heavy-light ground state mesons as a function of 1/Mg. 
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FIG. 24: Hyperfine splitting of the heavy-light-light baryons as a function of 1/Mg. 
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FIG. 25: M Ab — M§ as a function of 1/M B . 
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FIG. 26: — a4(A;,) + /J^(B) measured from the matrix elements is compared with the indirect 
measurement from the slope of mass difference M\ b — M#, which gives —0.03 ± 0.15 GeV 2 . 
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